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Gravitational lensing introduces the possibility of multiple (macroscopic) paths from an astrophys- 
ical neutrino source to a detector. Such a multiplicity of paths can allow for quantum mechanical 
interference to take place that is qualitatively different to neutrino oscillations in flat space. After an 
illustrative example clarifying some under-appreciated subtleties of the phase calculation, we derive 
the form of the quantum mechanical phase for a neutrino mass eigenstate propagating non-radially 
through a Schwarzschild metric. We subsequently determine the form of the interference pattern 
seen at a detector. We show that the neutrino signal from a supernova could exhibit the interference 
effects we discuss were it lensed by an object in a suitable mass range. We finally conclude, however, 
that - given current neutrino detector technology - the probability of such lensing occurring for a 
(neutrino-detectable) supernova is tiny in the immediate future. 
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1. INTRODUCTION 

Spacetime curvature allows, in general, for there to be 
more than one macroscopic path from a particle source to 
a detector. This means that there is a quantum mechani- 
cal interference phenomenon that may occur - at least in 
principle - with gravitationally- lensed, astrophysical neu- 
trinos that is qualitatively different from 'conventional' 
neutrino oscillation. The possibility for this different type 
of interference arises because - with, generically, each 
path from source to detector having a different length 
- a phase difference may develop at the detector due to 
affinc path difference(s). This is to be contrasted with flat 
spacetime neutrino oscillations which arise because differ- 
ent mass eigenstates generically have different phase ve- 
locities. One might expect, in fact, that gravitationally- 
induced neutrino interference ('GINf') exhibit a phe- 
nomenology partially analogous to that produced by a 
Young's double slit experiment, viz., regular patterns of 
maxima and minima across a detected energy spectrum. 
As we show below, for ultra-relativistic neutrinos, each 
maximum and minimum at some particular energy is 
characterised by, respectively, enhancement and deple- 
tion of all neutrino species (not relative depletion of one 
species with respect to another which characterizes flat 
space neutrino oscillations). 
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Below we shall provide the theoretical underpinning to 
all the contentions made above. We also sketch a proof- 
of-principle that this interference effect could actually be 
seen in the neutrinos detected from a supernova given 
a suitable lens. There are other situations where the 
GINI effect might, in principle, also be evident. Reluc- 
tantly, however, we conclude that pragmatic considera- 
tions mean that GINI effects will be very difficult to see 
in these cases. 



2. SURVEY 

Particle interferometry experiments enjoy a venerable 
lineage and - apart from their intrinsic interest - have 
often found utility in the measurement of intrinsically 
small quantities. The idea that the effects of gravity - 
the epitome of weakness as far as particle physics is con- 
cerned - on the phase of particles might become manifest 
in interferometry dates to the seminal, theoretical work 
of Overhauser and Colella 0. It was these researchers 
themselves, together with Werner 0, who were the first 
to experimentally confirm the effect they were predicting 
(in what has come to be labeled a COW experiment after 
the initials of these researchers: see Ref.Q for a review). 

Another interesting idea involving gravitational effects 
on interferometry of neutral particles - though, to the au- 
thors' knowledge, without yet having received experimen- 
tal confirmation - is the idea that gravitational micro- 
lensing of light might realise a de facto Young's double 
slit arrangement. There is an extensive literature devoted 
to this idea (see Refs. |llllillllli|li|ll|l2), which 
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has been labeled 'femtolensing' because of the natural 
angular scales involved for cosmologically-distant sources 
and lenses Femto-lensing is somewhat more closely 
analogous to the idea we present (indeed, as we show 
below, the analogy becomes exact in the massless neu- 
trino limit) than COW-type experiments. This is be- 
cause in femtolensing gravity not only affects the phase of 
the propagating photons, but is also itself responsible for 
the 'bending' of these particles so that diverging particle 
beams (or, more precisely, wave packets) can be brought 
back together to interfere. Furthermore, while the inter- 
fering particles are relativistic in both the femtolensing 
and GINI cases, they are non-relativistic in COW exper- 
iments. 

As far as sources go, light from GRBs has receivedpar- 
ticular attention in the context of femtolensing 0, S 13 • 
While we would, of course, also require astrophysical ob- 
jects as the sources for a GINI 'experiment', the sources 
best able to offer a chance for the detection of this effect 
are probably supernovae. A Galactic supernova would 
generate excellent statistics in existing solar (and other) 
neutrino detectors (thousands of events ~ see below). 
And with a much larger, generation of neutrino detec- 
tors on the drawing board - some having as one of their 
chief design goals the detection of neutrinos from super- 
novae occuring almost anywhere in our Local Group - 
prospects for the detection of the effect we predict can 
only improve with time [74] . 

By way of a pedagogical detour, please note the fol- 
lowing: we believe the 'time-delay' nomenclature is mis- 
leading in the context of either femtolensing or GINI ef- 
fects. It is much better, we contend, to think in terms 
of path difference(s). The idea of a well-defined time- 
delay belongs to classical physics. The time delay is - in 
the frame of some observer - the time elapsed between 
the arrival of two signals. These should have their ori- 
gins in the 'same' (macroscopic) event at a source, but 
then travel down different classical geodesies from source 
to detector. Now, from the viewpoint of quantum me- 
chanics, there is a limit in which the classical description 
just given makes sense and is useful. This limit is that in 
which the size of the wavepackets describing the signaling 
particles is small in comparison to the affine path length 
differences between the different classical trajectories un- 
der consideration. This limit will usually be satisfied in 
observationally-interesting cases of gravitational lousing. 
But this limit must not be satisfied if femtolensing or 
GINI effects are to be observed. Indeed, we require the 
opposite situation to pertain, viz, an affine path length 
difference of the order or smaller that the wavepacket 
size. This is required so that wave packets created in the 
same (microscopic) event can overlap at the detector - 
with interference effects being the result. In this sense, 
there is no time delay because the wavepackets have to 
be overlapping at the detector position at the same (ob- 
server) time, i.e., overlap must be satisfied at the space- 
time location of the observation event. Note further that, 
to paraphrase Dirac, each photon - or neutrino - only 



interferes with itself. So it is the wavepacket of the sin- 
gle particle that results from a single (microscopic) event 

- like the decay of an unstable parent particle - that, 
in simple terms, splits to travel down all the classical 
geodesies from source to detector, only to interfere when 
recombined there. The idea we are describing, then, does 
not require some weird (and impossible) analog of a 'neu- 
trino laser'; it works at the level of individual, particle 
wavepackets. 

Another strand that will be peripherally drawn into 
this paper is the behavior - at a classical level - of neu- 
trinos in a curved spacetime background (i.e., gravita- 
tional lensing of neutrinos treated as ultra-relativistic, 
classical particles). This topic became of immediate in- 
terest with the detection of neutrinos from SN 1987A 
|l3l H^ . Timing information from the nearly simulta- 
neous detection of these neutrinos and the supernova's 
photon signal and from the time and energy spread of 
the neutrino burst alone has been investigated in many 
papers as an enipirical limit on the neutrino mass scale 
(see, e.g., Refs. [Hill and Refs. [Ulilpj for reviews 
and the seminal references concerning this idea) and also 
as a probe of the equivalence principle over intergalactic 
distance scales |2y, |21| [73 . More speculatively but ger- 
mane to this work, the apparently bi-modal distribution 
of SN 1987A neutrinos observed by the then-operating 
Kamiokande solar neutrino observatory was given an ex- 
planation in terms of an intervening gravitational lens (in 
the 5 X 10^ Mq mass range: "2^). Further, the idea that 
astrophysical neutrino 'beams' might be gravitationally 
focused by massive objects like the Sun has been inves- 
tigated and it has been found that such focusing can 
amplify an intrinsic neutrino signal by many orders of 
magnitude (see ^| and 01 for more recent work). 

An early and important work treating the quantum 
mechanical aspects of neutrino propagation through a 
curved metric is that of Brill and Wheeler Their 
work is particularly important for its elucidation of the 
formalism that allows one to treat (massless) spinor fields 
under the influence of gravitational effects (i.e., the ex- 
tension of the Dirac equation to curved spacetime) 

As presaged above, in this paper we shall be particu- 
larly concerned with the phase of neutrinos (more partic- 
ularly, neutrino mass eigenstates) in curved spacetime. 
The seminal work treating the phase of quantum me- 
chanical particles in curved spacetime is that of Stodol- 
sky I2I. In this work the author argued that the phase 
of a spinless particle in an arbitrary metric is identical 
with the particle's classical action (divided by h). Later 
work conducted on neutrino oscillations in curved space- 
time mil m, has often - though not always Ha, 113 

- implicitly assumed the correctness of Stodolsky's con- 
tention (that the phase is given by the classical action) for 
spinor fields as well and taken this as its starting point. 
Somewhat ironically - as we set out in detail below - 
recent researches [sj have revealed that the equality of 
classical action and phase holds for spin half particles, 
but not for spin zero or one particles, or, at least not 
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in an unqualified sense. In any case, that Stodolsky's 
contention holds for spinors means a considerable sim- 
plification for our calculations as we can avoid directly 
treating the covariant Dirac equation. 

A full review of the literature (see [HIH IH 113, IH 
|3l|3l|3l|3l|3l|3l|3l|3i)on neutrino phase in the 
presence of gravity is beyond the scope of this work. Suf- 
fice it to say that most work here to date has been con- 
cerned with the calculation of neutrino phase in radial 
propagation of neutrinos through stationary, spherically 
symmetric, spacetimes. There is active controversy in 
this context as to at what order in neutrino mass (m^ 
or mf,?) gravitational corrections show up in the phase 
1^ ISa. l3^ . The answer to this hangs critically on 
how energy and distance, in particular, are defined [33| . 
We shall have to treat such issues carefully, but all the 
subtleties of this debate need not particularly concern 
us. This is because we are primarily interested in gravity 
not for its effect on phase per se, but for its ability to 
generate multiple macroscopic paths from source to de- 
tector. And it is what might actually be measured at the 
detector that concerns us. Detectors count the neutrinos 
- registered in terms of fiavor and (local) energy - that 
interact within their volume. From these one can infer 
interference patterns, but one does not, of course, have 
any direct experimental access to the phase difference(s) 
(a point that does sometimes seem to be forgotten). 

In regard to interference phenomenology, note the fol- 
lowing: whereas interference patterns with flat space neu- 
trino oscillations take the form of variations in neutrino 
fiavor ratios across energy, with GINI, because there will 
be constructive and destructive interference between the 
multiple allowed routes, one (also) expects to see, in gen- 
eral, maxima and minima (distributed across energy) in 
the counts of all neutrino flavors. These maxima and 
minima will be present irrespective of what measure of 
distance, say, we settle on (though, of course, they may be 
undetectably small in amplitude - but that is a separate 
issue). To put this in a different way, fiat space neutrino 
oscillations modify the relative abundances of neutrinos 
expressed as a function of energy whereas GINI effects 
can modify absolute abundances. 

Interestingly, of all the papers devoted to the topic of 
gravititationally-affected neutrino phase, to the authors' 
knowledge, only one |23| has previously examined GINI, 
which, to reiterate, is the idea that neutrino oscillations 
in the presence of gravitational lensing - or, to be strict, 
gravitational focusing - might present interesting phe- 
nomenology. (This is the analog of the femto-lensing de- 
scribed above that involved light.) To examine this idea, 
the authors of Ref. were obliged to develop a for- 
malism to deal with non-radial propagation of neutrinos 
around a lensing mass, and we shall adopt much of this 
formalism in the current work. Unfortunately, Ref. 0| 
contains an incorrect result which it is one of the major 
aims of this paper to point out. Moreover, other works 
which have considered gravitationally-affected neutrino 
phase contain results - and commentary thereon - which. 



if not strictly incorrect, can be misleading if one does not 
realise the restricted nature of their tenability. In brief, 
most authors have failed to consider the possibility of 
multiple paths. Any result which suggests the vanishing 
of neutrino phase difference in the massless limit [see, 
e.g., Eq. (13) of Ref. [H or Eq. (4 J] of [p] should be 
interpretted with extreme caution |76|. 

Essentially, the incorrect result in Ref. [23|, as al- 
luded to above is, then, that the phase for a neutrino 
mass eigenstate, j, propagating non-radially through a 
Schwarzschild metric is purely proportional to its mass, 
ruj, squared [see Eq. (58) of |27j and also Eq. (25) of ^]. 
This result is incorrect |77j: in the massless limit, the neu- 
trino phase in curved spacetime should reduce (modulo 
spin-dependent corrections which vanish for radial trajec- 
tories ^25] and are negligible except in extreme, gravita- 
tional environments 32]) to the result for photons. And 
the photon phase is not zero (otherwise the interference 
fringes - in space or energy - predicted by the femto- 
lensing literature would not be produced), even though, 
of course, the classical action is zero along null geodesies. 
Indeed, the photon phase is essentially proportional to 
energy. As we show below, furthermore, this oc E term 
is the leading order term for the neutrino phase as well. 

What has gone wrong when one's analysis misses the 
cx E term in the neutrino phase is that one has tried 
to simultaneously employ two incompatible notions: the 
fundamentally wave or quantum mechanical idea of phase 
with the particle notion of trajectory so that x is given in 
terms of t or vice versa. Even in the simpler case of fiat 
space oscillations, the introduction (often implicitly) of 
the idea of a trajectory - or, more particularly, a group 
velocity - into the calculation of neutrino phase leads 
to error (in particular, the recurring bugbear that the 
conventional formula for the neutrino oscillation length 
is wrong by a factor of two: see calcula- 
tions set out below, we show the reader how the error 
of introducing a trajectory can be avoided. Further- 
more, our method allows calculations to be performed 
along the actual (classical) paths [not trajectories; i.e., 
we have r(0), say, rather than r(t),(l>(t)] of the neutrino 
mass eigenstates, rather than taking the approach of cal- 
culation along the null geodesic employed in Ref ■ On 
the other hand, our approach also circumvents the obli- 
gation to introduce extra phase shifts 'by hand'. This 
artificial device becomes necessary when one offsets ei- 
ther the emission times or positions of the different mass 
eigenstates with respect to each other so that they arrive 
at the same spacetime point (see [28ij for an example of 
this). 

The plan of this paper is the following: in f^jwe de- 
scribe, for illustrative purposes, interference of neutrino 
plane waves propagating through flat space along both 
different (classical) paths and having, in general, different 
phase velocities. Then in f^Jwe describe the calculation 
of the neutrino mass eigenstate phase in a Schwarzschild 
metric, correcting an erroneous result that has existed in 
the literature for some time. We then set out, in ^Slthe 
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calculation of the analog of the survival and oscillation 
probabilities in flat space for neutrinos that have been 
gravitationally lensed by a point mass. In we exam- 
ine - at an heuristic level - questions of coherence that 
can effect the visibility of the GINI effect for neutrinos 
from supernovae. We give a proof-of-principle that the 
effect should be detectable. In fJHlwe describe some lim- 
itations of our method - which stem particularly from 
the assumption of exclusively classical paths - and set 
out improvements to be made in further work. Finally, 
in an appendix, we set out a wave packet treatment of the 
neutrino beam splitter toy model treated in §3 in terms 
of plane waves. We derive results here pertaining to the 
analog of the coherence length in conventional neutrino 
oscillations. 



NEUTRINO BEAM SPLITTER 



By way of an illustrative introduction to this topic we 
consider a toy model of interference effects that can arise 
when there are both multiple paths from a source to a 
detector, a la Young's double slit experiment, and dif- 
ferent phase velocities for the propagating particles, a 
la neutrino oscillations. Of course, interference requires 
that our experimental apparatus be unable to distinguish 
between the propagating particles (just as we likewise 
require for interference that the apparatus is unable to 
identify which path any single particle has propagated 
down). We therefore require the propagating and de- 
tected particles to be different objects. In this con- 
text, let us consider, for the sake of definiteness and 
relevance, a Gedanken Experiment involving an imagi- 
nary (flat space) neutrino beam splitter in the geom- 
etry illustrated in Fig. ^ We take it that the neutri- 
nos' paths can be approximated as two straight-line seg- 
ments along which momentum is constant in magnitude. 
We need only treat, therefore, one spatial dimension: 
J Pj.dx = J pjdx — Pj\x\^ where |a;| is the total distance 
along the two line segments. We expect that the quali- 
tative behavior of this device shall illustrate many of the 
features expected to emerge from interference of gravita- 
tionally lensed neutrinos. Note that for reasons of clarity 
we present only a plane wave treatment here, leaving a 
full wavepacket calculation for an appendix. We stress, 
however, that wave packet considerations - which allow, 
in particular, for a proper treatment of decoherence ef- 
fects - are, in general, important and must certainly be 
considered when one is dealing with neutrinos that have 
propagated over long distances (see §6). 

Let us write the ket associated with the neutrino 
flavour eigenstate a that has (in a loose sense) propa- 
gated from the source spacetime position A = {xA,tA) 



to detection position i? = (a;^, i^) as 
p j 

xexphi<i>^(i?,;L^^,T^^)]|^.,), (1) 

where = Ib — tA, L^^ is the distance from source 
position xa to detector position xb along one of a fi- 
nite number of paths labeled by p, ^^{Ej; L^^ ,T^^) = 
EjT^^ - Pj{Ej)Lp^ , with Pj{Ej) denoting the momen- 
tum of mass eigenstate j with energy Ej, and, finally, 
t/ is a unitary matrix relating the neutrino flavor eigen- 
statcs to the neutrino mass eigenstates. We have included 
the - ~ 

we should allow for a path-dependence to the amplitude. 

7^ Iq is the 



factor to account for the fact that, in general, 
d allow for a path-dependence to tl 
A situation where, for paths p and q, Ip 
analog of a Young's slit experiment in which the slits 
are not equally illuminated (thus reducing the visibility 
- see Eq. H92|l - of the resulting interference fringes). 
The Schwarzschild lens scenario explored in §§4 and 5 
presents a situation analogous to this: light or neutrinos 
propagating down the two classical paths from source to 
observer will not be, in general, equally magnified. Note 
that we choose throughout this paper - unless otherwise 
indicated - to work in units such that h = c = I G. 
The mass eigenstates are assumed to be on mass shell: 
Pj{Ej) = {Ej —rrVj)^^^ . The amplitude for a neutrino cre- 
ated as type a at the spacetime position A to be detected 
as type /3 at the spacetime position, B, of the detection 
event is then: 



J 



0j 



xexp[-i<i>^(i?,;L^^,T^^)]. (2) 

Assuming a stationary source, we can get rid of the un- 
wanted dependence on time by averaging over T^-^ in the 
above to determine a time-averaged oscillation probabil- 
ity analog at the detector position xb [4TL ■42i . ■43j . This 
maneuver gives us that 

\{iyi3Wa]XA,XB)\'^ oc / dT\{i^p\iya; A, B)\'^ 



oc 



pq jk 

xexp[i(p,(i?)L^^-pfc(i?)i^^)], (3) 

where we have Ej ^ Ef. ^ E because of the S{Ej — Ek) 
that arises from the integration over time. 

We find, then, after a simple calculation that the os- 
cillation probability analog becomes 

\{Uf3\iya;XA,XB)\'^ 

- ^ V ^P^Q E ^"i%^/3fe^a/c 



X exp 



, (4) 
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source, A 




FIG. 1: The geometry of either a double sht interference experiment or a two-image gravitational lens. The bold lines are 
the two neutrino paths, + and — , between the source, A, and the observer, B, and the geometry is defined both in terms 
of physical variables {va, tb, b± and s) and astronomical/lensing variables {Dd, Ds, Dds, 0± and /3). Also included is the 
deflection angle, a±{b), parameterised as a function of the physical impact parameter. It is implicit in this diagram that the 
deflector is rotationally symmetric, and that the paths are thus confined to the plane of the page defined by the source, deflector 
and observer. 



where the phase difference is given by 
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L 



2E 



AB 



(5) 



Note that in Eq. the normalization has been deter- 
mined by requiring that \{i'p\^'a]XA,XB)\^ < 1- The 
presence of the < sign (as opposed to a simple equality) 
comes about as the interference between states propagat- 
ing down different paths can result in minima at which 
the total neutrino detection probability is zero, in which 
case the usual unity normalisation is impossible (so now 
Kvfslva', xa, xb)\'^ no longer has a direct interpretation as 
a probability) . This is qualitatively different behavior to 
that seen in neutrino oscillations where, given maximal 
mixing between i/a and Vf3, an experiment might be con- 
ducted in a position where only neutrinos of type a, say, 
are to be found or, alternatively, only of type /J, but never 
in a position where none can be found in principle. 

On the other hand the behavior explained above - 
involving interference minima and maxima - is obvi- 
ously analogous to what one would expect in a double 
slit experiment or similar. In fact, we shall show be- 
low that the phenomenology of neutrino interference - 
when there is more than one path from source to detec- 
tor - is a convolution of the two types of interference 
outlined above. This means that, in general, one can- 
not simply re-cast \{i'p\i^a',XA,XB)\'^ in terms of a con- 
ditional probability, separating out the overall interfer- 
ence pattern (with its nulls, etc.) from the conditional 
probability that any detected neutrino has certain prop- 



erties. To understand this point, imagine setting all mix- 
ing angles to zero so that the mass eigenstates and weak 
eigenstates are identical. The point now is that the in- 
terference patterns for the various detected weak/mass 
eigenstates are still different: the phase difference (which 
now stems purely from path difference) is dependent on 
the mass of the neutrino species involved. In other words, 
K^^/glfa! xa, xb)\'^ will not, in general, factorise into a con- 
ditional probability multiplied by an interference enve- 
lope because the putative interference envelope is differ- 
ent for different mass eigenstates. As we show below, 
however, in the ultra-relativistic limit, the cross term 
(that mixes different path indices with different mass 
eigenstate indices) always turns out to be small with re- 
spect to the other terms. In this limit, then, factorization 
is a good approximation. 

We now consider two particular, illustrative cases of 
the neutrino beamsplitter gedanken Experiment for which 
we calculate relevant phase differences and oscillation 
probability analogs. 



3.1. Two Path Neutrino Beam Splitter 

A particularly perspicuous example is given by the two- 
path example of the above equations. For this we specify 
a reference length, L, to which the two paths, of lengths 
L-i_ and are related by 



L±=L± 



AL 



(6) 
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This means that there are four phase difference types as 
labeled by the path indices p,q (z {+, — }, viz: 

^ 6m%L + AL/2 
i. A$++ = + — ^- — 

6m%L - AL/2 

»■ A<f7.- = 

A<i>V = -E (l -"^"^ ^ f\ AL ^ 

ok I 4^2 I 



2E 



(7) 



In this case, then, Eq. becomes 

\{yi3\Va]XA,XB)\^ 
1 



/+ + /_+ 



^ Uc^jUfj^UpkKk exp -i 



2£; 



X < /_!_ exp 



.Sm%AL\ 



AE I 



/_ exp — i- 



AE 



2\Jl+I- cos 



/ TO" 



AE^ 



(8) 



rigorous, if simplistic, results derived above with heutris- 
tic arguments to investigate a more realistic double slit 
neutrino interference experiment. The choice of a double 
slit experiment is particularly relevant not only because 
of its links with more familiar interference phenomena, 
but also because a point-mass gravitational lens admits 
two (significant) paths from source to observer. Thus 
the results obtained in this section should provide a use- 
ful guide to the qualitative behaviour of GINI in curved 
space time considered in 50] 

The relevant geometry is illustrated in Fig.^ showing 
the two classical paths (+ and — ) from source to detector, 
with the important introduction of a more physically mo- 
tivated set of lengths than the reference length, L, used 
previously. The entire experiment is taken to be planar, 
with the coordinate system origin on the line defined by 
the two slits. (This is an arbitrary decision at present, 
but will coincide with the position of the deflector in SQl) 
The slits are defined by their positions, 6+ and 6_ (which 
will correspond to image positions in ^4.51) ; the source po- 
sition is given by both its radial coordinate, r^, and its 
perpendicular offset, s; the observer position is defined 
by its radial coordinate, r^, alone. There are several 
other plausible ways in which this geometry could be de- 
fined, but all derived results become equivalent under the 
assumption that \h±\ <ti applied throughout. 

Note also that the observer and source are interchange- 
able. 

As discussed in fj^l all flat-space interference phenom- 
ena can be treated in terms of path lengths. The lengths 
of the two paths illustrated in Fig. ^ are 



This is an interesting result. It shows that the interfer- 
ence factorises into a conventional, flat space oscillation 
term and an interference 'envelope' in curly brackets (as 
might be expected from a Young's slit type experiment). 
If we now further particularize to the case where /+ — I~, 
this envelope term becomes 



AE 



cos 



E 1- 



AE^ 



2 cos 



El 



\ AL 






cos 







(9) 

which obviously reduces to the expected photon interfer- 
ence term c>c cos^ in the massless limit. Later (in fj^J 
we shall see that this sort of factorization property also 
arises, under a different set of assumptions, with gravi- 
tational lensing of astrophysical neutrinos. 



{rA + tb) 



L± = {rl + bl) + [r\ -s^ + (&± - s)'] ' (10) 



± 



,1/2 
-sh± 



where, in most cases, h- and h+ are of opposite sign, and 
the second line explicitly utilises the fact that ^ 
^{A,B}- The path difference is thus 



TA + TB 



(11) 



2rB 



^rATB 

and the average path length is 
L ~ (rA + tb) 

1 + -^ bl + bl- 



rA + TB 



■s {b. 



b' 



(12) 



Atatb 



2rB 



TA + TB 



-S{b++b-) 



3.3. Schwarzschild slit geometry 



3.2. Double slit experiment geometry 

Having considered a hypothetical neutrino plane wave 
beam-splitter f H3.1|l . it is now possible to combine the 



In a laboratory-based double slit experiment the two 
slit positions can be chosen arbitrarily, but in the case of 
gravitational lensing the impact parameters of the beams 
are determined by a combination deflector and source pa- 
rameters (SQ. Given the Schwarzschild metric around a 
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point-mass f H4.5|l . the assumption of small deflection an- 
gles implies that the source position and impact param- 
eters are related by [cf. Eq. l(7^ ] 



bi-bi = Abi 



s{b+-b^). 



(13) 



Applying this result in the more general context of the 
double slit experiment, Eq. (|13|l can be rewritten as 



AL I _ = Li 



^rATB 



(14) 



Similarly, the mass-length expression that appears in, 
e.g., Eq. Q can be simplified to 



mlL. 



{rA + tb) 



(15) 



8m% — 

2rArB 



{mpl-mlbl) 



Substituting these expressions in Eq. (jSJ) then gives 
the phase difference between mass eigenstates j and k 
travelling down paths -I- and — as 



^E{rA+rB) 



Abi 



2rArB 
rA 



2E 



TB {mpl 



{rA + rB) 
-mlbl) 



2E 



2rArB 



(16) 



This expression includes contributions from both differ- 
ent phase velocities and different path lengths, and can 
be understood further by considering the special cases in 
which (1) different mass neutrinos travel down the same 
path or (2) the same mass eigenstate travels down differ- 
ent paths. 

1. From Eq. (O for the general case of the phase differ- 
ence between the same mass eigenstate propagating 
along different paths we find 



E ^- 

2E 



{Lp — Lq 



(17) 



So that if we further particularize, as above, to 
paths passing through slits at 6+ and b- we find 
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TO' 
E 

2E 



{ra+n)Ab'\ 
'2rarb 



(18) 



Notice in the above the similarity to the phase dif- 
ference for an ordinary Young's double slit type ex- 
periment using photons, namely. 



A$7 = -Ev{Lp~Lg) 



(19) 



where v is the (phase) velocity of the interfering 
particle (which we assume to be relativistic). We 
shall see below that the analog of this phase - es- 
sentially proportional to energy x path difference 
- has been missed in the existing literature on neu- 
trino oscillations in curved space. This has led to an 
incomplete result suggesting that the phase differ- 
ence vanishes in the massless limit even when there 
is more than one path from source to detector. 

2. Again for the general case, for the phase differ- 
ence between different mass eigenstates propagat- 
ing along the same path (i.e., the analog of the 
usual phase difference encountered in neutrino os- 
cillation experiments), we find from Eq. |(SJ) 



2E 



where x g {p, q}. 



2E 



[La +Lb){l + 



2LaLb 



(20) 



It is worth keeping the above expressions in mind when 
considering the results for the neutrino phase difference 
in curved spacetime presented in H4.5I As will be seen, 
the results obtained in this more complex physical situa- 
tion are analogous to those derived above, e.g., compare 
Eq. 10 with Eq. ^ and Eq. ^ with Eq. 



4. THE PHASE OF A NEUTRINO MASS 
EIGENSTATE IN CURVED SPACETIME 

A neutrino beam splitter is in the realm of fantasy 
- except for the interesting case of gravitational lensing 
of neutrinos: a gravitational field can bring to a focus 
diverging neutrino beams, and therefore provide for mul- 
tiple (classical) paths from a source to a detector. In the 
remainder of this paper we explore whether any interest- 
ing, quantum mechanical interference effects can arise in 
this sort of situation. 

We shall be concerned below, therefore, with deriving 
an expression for the neutrino oscillation phase in curved 
spacetime, in particular a Schwarzschild metric (this pro- 
viding the simplest case in which gravitational lensing is 
possible). Here we shall follow the development laid out 
in |23| and ^30] but, importantly, we shall also employ 
the prescription set out in f|31that allows for the removal 
of time from consideration in the oscillation 'probabil- 
ity' by integration over T = — where A = {va, tyi) 
and B = (tr, tp) are the emission and detection events 
respectively |41j . Note that we are assuming the semi- 
classical limit in which gravity is not quantized and its 
effects can be described completely by a non-fiat metric, 
9 til' 7^ Vfj-i^- 

The procedure we follow is to start with the gener- 
alization of the equation for a mass eigenstate's phase 
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in flat spacetime to curved spacetime first arrived at by 
Stodolskyi2^: 



where 



Pi'' 



P 



ds 



(21) 



(22) 



is the canonicaUy conjugate momentum to the coordinate 
xf^. Actually, as pointed out by Alsing et al. in Ref. [s^ . 
Stodolsky's expression for the phase is missing, in gen- 
eral, small correction terms that arise from quantum me- 
chanical modifications to the classical action. These vary 
according to the spin of the particle under consideration. 
Completely fortuitously, the would-be correction terms 
are identically zero in the case of spin half particles in a 
static metric (whereas for particles with, e.g., spin zero 
or one they are non-zero) so the Stodolsky expression 
happens to be exact for the Schwarzschild metric and 
many other cases of interest. Note in passing that this 
restriction to a static metric means that this technology 
cannot - as it stands - treat, e.g., particle phases in a 
cosmological context. 

We now introduce the metric of the Schwarzschild 
spacetime. This may be written in radial co-ordinates, 
x^^ = (i, r, 1?, (fi), as 



where 



B{r)di^-B{ry^dr^-r^dd^-r^ sin^ Mtp^ , (23) 



B{r) 



1 



2GM\ 



(24) 



and G is the Newtonian constant and M is the mass of 
the source of the gravitational field, i.e., the lensing mass. 
Given the isotropy of the gravitational field the motion of 
the neutrino mass eigenstate will be confined to a plane 
which we take to be the equatorial one, i9 = 7r/2 and 
dd = 0. 

The relevant components of the canonical momentum, 
Eq. are, then ^J: 



(k) 
Pt 



mkB{r) 



dt 

ds ' 



p«=-m,S-i(r)J, 



and 



(k) 



,d(p 
ds 



(25) 



(26) 



(27) 



These are all inter-related through the mass-shell condi- 
tion il: 



ml = g^'pi^'pi''' 



1 / (^)^2 

-l^(pp))^-s(.)(pW)^-%l 



(28) 



Given that the components of the metric are indepen- 
dent of the coordinates t and ip, the momenta associated 
with these quantities, p[''' and p^^' shall be conserved 
along the classical geodesic traced out by i^k- We define 
these constants of motion as Ek = p\''' and Jk = —p^p' . 
These two are, respectively, the energy and angular mo- 
mentum seen by an observer at r ^ oo for the A:th mass 
eigenstate They are not identical with the energy 

and angular momentum that would be measured for I'k 
at some definite, finite position r. In general, however, 
one may relate these quantities using the transformation 
law that relates a local reference frame {x"} — {t, f, i), (p} 
to the frame {x^} = {t, r, ■d,(p}^: 



5^. = L'^.L^^rj^,, (29) 



where the L",'s are the coefficients of the transformation 
between the two bases: 



L\ = VM, = V\9rr\, 



L\ = VM: L% = ^\g^^\. (30) 

So we have, in particular, that the local energy is given 
by M- 

Et'\r) = \gu\-^'^Eu - B{r)-^'^Eu. (31) 



4.1. Calculating the Phase Difference 

Given the above definitions, we now have that: 



$^(S,A) 



[Ejdt - pj{r)drp - Jjdipp], 



(32) 



where we have implicitly defined Pj{r) = ~~pr . Note 
that we have explicitly introduced the path index p which 
allows for the possibility of multiple paths from source to 
detector. Again, however, the integration over t is inde- 
pendent of the path as the endpoints of this integration 
are defined by the emission event and detection events. 
In fact, as discussed above, Ej is conserved over classical 
paths, so that if mass eigenstate j is assumed to travel 
down such a path, we can calculate the phase it accumu- 
lated after leaving the source to be 



$^(i?,^) 



Ejdt 



dr„ 



Ej{tB-tA) 



Pj{r) + Jj 



dip 
dr 



drp. (33) 
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Of course, the quantity that governs the oscillation phe- 
nomenology is the phase difference A$^^ where, gener- 
ically, interference can be between different mass eigen- 
states and/or different paths (cf. discussion in As 
things stand this quantity would be parameterized in 
terms of both t and r: 



A^PlirB,tB,rA,tA) = {Ej ^ Ek){tB - tA) 

dr„ 



dr„ 



(34) 



We therefore follow the prescription set out in jlj to rid 
ourselves of the unwanted time parameter: we assume 
a stationary source and integrate the interference term, 



exp 



iA$ 



over the unknown emission time tA (or, 

equivalently, the transmission time T = Ib — tA)- This 
results in a very useful 5{Ej — Ek). 

Note here that though the energies of different mass 
eigenstates are different - so that the S {Ek — Ej) 
arising from the time integration would seem to imply 
no interference - in fact, in a correct treatment, massive 
neutrinos are described by wave packets, not plane waves 
as here. This means that, though the average energies of 
different mass eigenstate wave packets are, in general, 
different, each massive neutrino wave packet has an en- 
ergy spread and the detection process can pick up the 
same energy component for different massive neutrinos 
(see Refs. plll47j ). If the energy spread of the wave pack- 
ets is small there is a suppression factor that, formally, 
can only be calculated only with a wave packet treatment 
(cf. 33; but, can also be assessed at the heuristic level 
(cf. f©. 

Let us see how all the above works in practice. 



4.2. Radial Propagation 

We consider first the simple case of radial propagation, 
in which case there is a single classical path from source 
to detector. Along this path, the angular momentum 
vanishes and we have: 

A^jk{rB,tB,rA,tA) = (Ej - Ek){tB - tA) 



[pj{r) - pk{r)]Ar (35) 



We can determi ne P j(r) — Pkir) from the mass-shell re- 
lation, Eq. if^ 



me Pi 

m 

Pk{r)=±^^^[E^[^B{r)^,, 



(36) 



where the -I- sign refers to neutrinos propagating out- 
wards from the gravitational well and the — sign to neu- 
trinos propagating inwards. We can further simplify this 



relation by employing the binomial expansion which, as 
in the flat space case, holds for relativistic particles: 



^El-B{r)mi ^ E^ - B{r)^^, (37) 

where Eq is the energy at infinity for a neutrino mass 
eigenstate in the massless limit (see |27] for a detailed 
account of the region of applicability of Eq. (|37|l ). We 
therefore have that 

The phase difference then becomes 



I^^jk{rB,tB,rA,tA) 

- {^3 ~ ^k) 



{tB-tA)T / 

J TA 



dr 

W) 



Given the oscillation 'probability' shall be, following our 
previously-establsihed procedure (cf. Sj^J , integrated over 
T = tB — tA, the relevant phase difference can be seen to 
be 



A$jfe(rB,rA) 



2E, 
2Eo 



TB - rA\ 



(40) 



where E^ — Ej — Ek and Eq is the energy at infinity 
for a massless par ticle and, as in flat space, the following 
relation holds l27ll: 



Ek ^ Eo 



0{ ^ 

2Eo 



(41) 



To digress a little, note that the result presented in 
Eq. (gni), arrived at previously [Ulli, HE Ell , must be 
interpreted with some care: in Eq. H4()|l one must keep 
in mind that the radial distance jr^ — is a coordinate 
distance, and not the proper distance the various mass 
eigenstates experience (except in the flat space case to 
which Eq. (|4(J|I clearly reduces in the limit of a van- 
ishing lensing mass) and that Eq does not represent a 
locally-detected energy. Following j23|, however, we can 
convert the phase difference so that it appears in terms 
of these parameters. The proper distance is given by (cf. 
Eq. ^): 



^prop 



/.g^rdr 



rB-rA + GMln—, 



(42) 
(43) 



where in the second line we have assumed the weak field 
limit holds. This demonstrates that, in a gravitational 
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field, the length relevant to the calculation of phases, 
\rB — t'a\i is actually shorter than the distances experi- 
enced by the propagating particles, Lprop- Substituting 
Eqs. (EU and (O into (gOJ we determine that (cf. 0): 



A$jfe(rs,rA) 



Ar7i|j,Z/pi.op 



1 - GM 



4.3. Non-Radial Propagation 



(44) 



In the above we have also employed the fact that the 
angular momentum of mass eigenstate j (traveling along 
path p) at infinity is given in terms of Vj^s energy at 
infinity, Ej^ the impact parameter along the path being 



considered, 



bp, and Vj 



's velocity at infinity, v 



(oo) 



Jj =EjbpVj 



(oo) 



Ejbp 



2El 



(48) 



We turn now to the more interesting case presented by 
non-radial propagation. Here there will be, generically, 
more than one path for the mass eigenstates to take from 
source to detector and we have the possibility, therefore, 
of interference between particles on these different paths. 

The phase difference we must calculate is given by 
Eq. To proceed with this calculation we must de- 

termine a value for 



Pjir) 



We can further evaluate Eq. 07| by replacing the path- 
dependent impact parameter, bp, with the minimal radial 
co-ordinate for the same path, rg. The relation between 
these two is found by noting that at the position of closest 
approach the rate of change of the co-ordinate r with 
respect to the angle (f> vanishes [23| . For the massive case 
fEg. I45|) . this implies that: 



Firstly recall that Jj is constant along the classical path 
taken by Vj. Now, using the fact (see Eq. (101.5)) 
that 



I ™-j 2GM 



(49) 



dr 



±- 



-^1 



Bir) 



B{r) 



^P(i-f) 



(45) 



and, given that from the mass-shell relation we have [see 
Eq. (49) of 123] : 



Bir)p,{r) ^ ±dl-B{r)% 



B{r){l-biy) 



1 



B{r)bl/r^ 2E] 



(46) 



Employing Eq. (|49|l . taking the weak field limit, and also 
expanding to 0{ni?JE'^) we find that: 



dr 



E, 



r 



GM- 



2r - 3rg 



GM- 



Vr2-(rg)2(r-frg) 



(50) 



where bp is the impact parameter for path p, one may 
determine that: 



± 



E, 



B{r)^l-B{r) 



m| l + (l-2B{r))^ 
2E] 



(47) 



With this result in hand, we can complete the calculation 
of Eq. |(3S|), the phase accumulated by mass eigenstate 
in non-radial propagation from spacetime position A = 
ijA^tA) to B = {rB,tB), where either or tb is the 
minimal radial co-ordinate encountered over the journey 
(i.e. the path is either non-radially inwards or outwards 
but not both). After an elementary integration we find 
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that 



2GM In 











/ rs 









t2 r 



± 



2Ei 



P\2 



-GM 









-rl\ 


1 rs 









(51) 



where the upper signs pertain if dr is positive (outward 
propagation) and the lower if dr is negative. 



4.4. Neutrino Lensing 

Finally let us consider the case of gravitational lensing 
of neutrinos. In this case the neutrinos propagate non- 
radially along classical paths, labelled by index p, from 
radial position ta, inwards to a path-dependent minimal 
radial co-ordinate Tq, and outwards again to a detector 
situated at radial co-ordinate rs- As presaged above, in 
this situation there will be (at least potentially) interfer- 
ence not only between different mass eigenstates prop- 
agating down the same classical path, but also between 
mass eigenstates propagating down different paths {p and 
q, say). Taking into account the sign of the momentum 
along these two legs, we find, following the developments 
above, that the relevant phase is given by 

<i>?(B,^) 

^Ej{tB~tA) 



E, 



^Ai^+2GMinf!:A±V5ZEI!^ 



-GM 



TA + rl 



+ {ta rs) 



2E, 



GM 



(52) 



where (r^ —^rg) mean add another term of the same 
form but with ta replaced with rs- 



Before proceeding any further with the calculation it 
behooves us here to establish the plausibility of Eq. 1)52(1 
by showing its relation to results known from some sim- 
pler cases. In the M ^ limit this equation becomes 




Now refer back to Fig. ^ and take the coordinate ori- 
gin on the diagram to denote the position of a lensing 
point mass. In the massless case, the two classical paths 
reduce to the single 'undeflected' path denoted by the 
dashed line in the diagram. Denote the minimal ra- 
dial coordinate along this path by vq (which intersects 
the dashed line at right angles). Clearly, then, the ge- 
o metrical leng th of th e path from source to detector is 
y/r\ — Vq + \/t\ — Tq . Now, given we know that in flat 
space the phase of mass eigenstate fc is given by Eq. 1)21(1 
with Minkowski metric, viz: 



= 

then, for the M 
comes 



Ej{tB - tA) - Pj-ixB - xa). (54) 
case illustrated in Fig. this be- 



E,{tB-tA)-E, 1--^ 




(55) 

With Eq. H55() we have, then independently established 
the plausibility of Eq. ((53() , once one takes into account 
the fact that, in the massless lens case, all classical paths 
converge on the same undeflected path (as mentioned 
above) so that in this limit Tq = rg = ro Vp, q. 

The other limit of interest is to take rrij — > in 
Eq. 1(52(1 . In doing this - and then setting the temporal 
and spatial contributions to the phase equal as appropri- 
ate for a null geodesic - we find that we have re-derived 
the Shapiro time delay [see, e.g., Eq. (8.7.4) of Ref . [s^]. 

Continuing with our main calculation, we can re-write 
Eq. I(52() in terms of bp by inverting Eq. ((49(1 . If we also 
expand to 0{bp/r'^^j^), we find that 

— Ej{tB — tA) 

-Ej{rA + rB) 



2GM 



^rATB TA + TB 



rB 



2GM \ 



^taTb ta + rB 



^1 



(56) 
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We can now find the phase difference, which allows for in- 
terference between different paths and/or different mass 
eigenstates, by the usual integration over T (so that we 
have El, ~ Ej = Sfe): 



^ +E^{rA 



2E„ 




where Afo^^ 



(57) 



A$^^(rs) 



and, in our notation 
denotes the phase difference between mass eigenstate ] 
traveling down path p and mass eigenstate k traveling 
down path q |78l |. Eq. I)57|l is one of the major results 
of this paper. Note that the presence of the cx Ey term 
in this equation - missed in Ref. [27| - ensures that the 
phase difference behaves properly in the massless limit 
(i.e., does not vanish). In passing, also note that the 
above equation satisfies the discrete symmetry of swap- 
ping B and as it should: the same result must be 
obtained for the phase difference (in a static spacetime) 
if we swap the positions of source and observer. 

Also recall that, excluding the case of perfect align- 
ment of source, lens, and observer, there are only two 
possible classical paths from source to observer for the 
Schwarzschild case. These we label by -I- (this path hav- 
ing an impact parameter somewhat greater than the im- 
pact parameter for an undeflected ray) and — (this path 
having an impact parameter on the 'opposite' side of the 
lens to the undeflected ray). We require, therefore, that 
p,q & {+, — } and in the particular case that we are con- 
sidering interference between the same mass eigenstates 
propagating down different paths Eq. (|57ll becomes 



A$+- {rB,rA) 



+E^{rA + rB)i + 



AGM 



2rArB TA + TB 



m] Abl_ 
2E„ 2rArB 



In 



6_ 



(58) 



Alternatively, in the case of different mass eigenstates 
traveling down the same path (i.e., 'ordinary' neutrino 
oscillations, but in curved space), Eq. lf57|l becomes 



A<^'U^B,rA) 



2E, 



{rA + tb) 



^tatb 



2GM 

TA + TB 



(59) 



Note that Eqs. and give us that 

A^f^irs^TA) = A$PJ(rB, r^) + Aci>^^ (rB,rA) 

ir.^rB)^ 
2rArB 



O 



2E, 



(60) 



This correction term will be small with respect to other 
terms (given our assumptions of ultra-relativistic neutri- 
nos and undeflected impact parameters small with re- 
spect to the overall distances between source-lens and 
lens-observer). In fact, the third term of Eq. H57|) can be 
expected to be suppressed with respect to the first term 
by 0{m'^ /El) and with respect to the second term by 
O {iP' / {r AT b)) ■ The consequence of this is that the phase 
may be written 



satisfying what we label 'separability', where 



(61) 



A^P'i = E,{rA+rB, 



pq 



AGM 



\2rArB ta + tb 



In 



2e:^'^ 



tb) 



pq 



2rArB 



(62) 



and 



A$ 



jk 



'ok 



2E, 



(rA + TB) 1 



2GM 



2rArB rA + tb 



where 



— Y 



2 j and b = 
^ N, 



path 



(63) 



(64) 



with Ni, the number of neutrino mass eigenstates and 
A^path the number of classical paths from source to detec- 
tor (two in the case of the Schwarzschild metric). What 
Eq. H61|l says in words is that the phase difference that 
develops between source and detector is due to two effects 
that can be considered separately: (i) a phase difference 
~ independent of which mass eigenstate is under consid- 
eration - that develops because of the different lengths of 
the paths involved and (ii) the phase difference that de- 
velops because the different mass eigenstates travel with 
different phase velocites. This situation is analogous to 
two runners who run along two very similar - though not 
identical - paths, with similar - though not identical - 
velocities: to first order, the difference in the finishing 
times between the two depends on terms proportional 
to the difference in lengths of the two courses, AL, and 
the difference in the runners' velocities. Aw, but not, by 
definition, on terms oc AL Aw. 
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4.5. Phase Difference in Terms of Conventional 
Lensing Parameters 

To facilitate interpretation of the above results in an 
astrophysical context - and, eventually, to introduce an 
evolving cosmological model - it is useful to re-express 
the phase difference in the language of standard gravita- 
tional lensing theory (despite the fact that the particles 
being lensed are not photons). 



4-5.1. The lens equation 

The classically allowed neutrino paths in the presence 
of a deflector can be derived by reconsidering the ge- 
ometry shown in Fig. ^ Under the assumption that 
\b±\ ^ r^A.B], the source offset, s, can be related to 
the impact parameter, 6, by the lens equation: 



-b + TA a{b), 



(65) 



where a{b) is the deflection angle of the lens as a function 
of impact parameter. 

It is standard practice to reexpress the lens equation 
in terms of angular variables. This entails replacing the 
source offset and impact parameters with angles (on the 
sky of the observer) and radial coordinates with line- 
of-sight distances. These conversions are summarised 
graphically in Fig.Q] which leads to the following replace- 
ments: rA — > {D'^s + s^Y^'^ ~ Dds, where Dds is the dis- 
tance from deflector to source; rs — > Dd, where Dd is the 
distance from observer to deflector; and ta + tb — > Dg^ 
where Ds is the distance from observer to source. The 
notation employed for the distance measures is sugges- 
tive of their being the angular diameter distances used 
to relate angles and lengths in an evolving cosmological 
model, and they fulfil an analogous role here. It is most 
important to note, however, that they are not true angu- 
lar diameter distances and the following results are only 
quantitatively valid on scales sufficiently small that the 
expansion of the Universe can be ignored (e.g., the Milky 
Way or the Local Group) . These results will be extended 
to an evolving cosmology in |54| . 

The above caveats notwithstanding, the angular posi- 
tion (relative to the deflector) of an image with impact 
parameter b is now simply 



b 

TB 



b_ 



(66) 



and the position of the source can be given in terms of 
an unobservable angular parameter (3 as 



TA + TB 



S 



(67) 



Inserting these definitions into Eq. 1)65(1 , the lens equation 
becomes 



/3: 



Dd 



-aiDdS 



(68) 



The position(s) of the images formed by a source in a 
given position can then be found for a given choice of 
deflector model. 



4.5.2. The Schwarzschild lens 

In a Schwarzshild metric, the total angular deflection 
of a particle of mass m impinging on a point-mass M 
with undefiected impact parameter b is (see, e.g., [s^): 



n(6) 



AGM 1 
b~2 



(69) 



where v^o is the particle's speed at an infinite distance 
from the mass and it has been assumed that b ^ 2GM — 
Rs, the defiector's Schwarzschild radius. For an ultra- 
relativistic particle, this becomes 



areiib) 



AGM 



(70) 



where E is its coordinate energy (equal to the energy 
measured at an inflnte distance from the mass). 

For astrophysical neutrinos, however, m? /{2E'^) ^ 1 
and so it is an excellent approximation to assume they 
travel along classical photon paths, for which 



aiight(&) 



4GM 



(71) 



Previously we have been rigorous in taking the classi- 
cal paths of massive particles from source to observer 
but, as will be seen below, this assumption is entirely 
self-consistent when dealing with weak-field gravitational 
effects and ultra-relativistic particles. Note also that 
a corollary of this approximation is that the different 
mass eigenstates are assumed to travel down identical 
paths (whereas in reality the heavier eigenstates will fall 
marginally deeper into the deflector's potential well). 

Applying the above deflection law to Eq. gives the 
point-mass lens equation as 



(3^9--^ 



where 



= ,/4GM 



Dds 
DdD, 



(72) 



(73) 



is the Einstein radius of the lens. This is the angular 
radius of the circular image that would be formed in the 
case of perfect source-deflector-observer alignment (i.e., 
f3 — 0) and thus depends on distance factors as well as 
the lens mass. Solving the lens equation then gives the 
image positions as 



4^1 



(74) 
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This also implies the useful Schwarzschild-specific result 
that 



Mi 



(75) 



Having found a relationship between the angular posi- 
tion of a neutrino source and its images, the expression 
for the phase difference given in Eq. H57II can be recast in 
a form containing only line-of-sight distances and angu- 
lar variables. This yields the Schwarzschild-specific result 
that 



ds 



3l In 



5m 



DAl- 



1 DriD.^m^ 



2E, Dds 



(76) 



The second term in this equation is simply the phase dif- 
ference that develops between mass eigenstates j and k 
traveling along the same path for distance Dg, with a 
small correction for the presence of the deflector. The 
first term encodes the path difference along the trajec- 
tories p and q, with separate contributions from the ge- 
ometrical effect (oc A0p^) and the reduced coordinate 
velocity close to the deflector [cx ln(|0p/0g|)]. The final 
cross term is the leading order contribution from differ- 
ent eigenstates traveling down different paths. From the 
discussion in the previous section, this term will be small 
in general. 

Given that interference effects can only ever be im- 
portant when the detector cannot resolve different image 
positions (i.e., it cannot know down which path the neu- 
trino has travelled), having the phase difference in terms 
of 9p and 9q is not as useful as expressing it as a function 
of the (angular) source position, (3. 

For the Schwarzschild lens the conversion from 9 to f] 
is given in Eq. H72|) . and substituting this into Eq. I|76l) 
then gives (for mass eigenstate j down path + and mass 
eigenstate k down path — ) 



E, 



DdDs 



A91 



2E, 

1 DdDs I 
2E„ D,s 4 



9|ln 



fi - v//32 + ^91 



5m%{(i^+29l) 



(77) 



Thus the phase difference is expressed in terms of ess- 
ntially independent astronomical variables: the line-of- 
sight distances between observer, deflector and source, 
the mass of the deflector (encoded uniquely in 9-^ once 
the distances have been chosen) and the perpendicular 
source offset, (3. 

Most of the important results obtained towards the end 
of ^4.41 can be recast similarly in terms of standard lens- 
ing variables, either in terms of the unobservable image 
positions or the source position. Assuming separability 
(see ii4.4l) . for instance, the contribution to the A<I> due 
to path difference effects alone (Eq. I62|) can be written 
as 



which for p 



E, 



DdDs 
D,s 
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DdDs Ml 
2Ei, Dds 

- and q = 



2 

becomes 



(78) 
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(79) 



Similarly, the contribution due solely to the different 
phase velocities of two mass eigenstates traveling down 
the same path (Eq. 1^ can be expressed as 



Sr 



2E,, 



Ds 



(80) 



5. THE OSCILLATION 'PROBABILITY' 

With the above results, we can now calculate the ana- 
log, in curved space, of the flat-space neutrino oscillation 
probability: 



(81) 



\{vp\v^,Ds)\'' (X j dT\{vp\v^-A,B))\'' 
so that we can write 
\{vp\v^,Ds)\'' 

pq jk 

X exp [-iA$P^(i?,) 



(82) 



where Ip and Iq account for the fact that different paths 
may be differentially magnified by a lens. We remind the 
reader that | {vp\va, -D^) ^ is no longer strictly a probabil- 
ity - see ^ 
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In the case that the 'separability' defined by Eq. H61() 
is satisfied, if, for the moment, we are interested only 
in determining the (energy) spacing of the interference 
maxima and minima, we need only consider a plane- 
wave-like calculation (and can therefore set to one side 
the coherence length effects and so on that emerge from 
a wavepacket calculation). So, following considerations 
similar to those that lead to Eq. H124|l in the appendix 
we can calculate that 

\{vf,\vo..Ds)\'' 



iVpath/ + 2 VV«cos(A$P'') 



p.q<p 



sm ' — 



V 2 



-sin(A$jfe) 



(83) 



where A'^path is again the number of classical paths from 
source to detector (two in the case of a Schwarzschild 
metric), / = ^p/-^path, and the normalization, \N\'^ is 
again given (cf. ^ by requiring that J2f3 IWpWa, Dg) P < 
1, i.e., 



Eq. H83() then becomes: 



(3^ + 40\ 



■{/32 + 2(?|[l + 2cos(A$+~)]} 



j,k<j 



jk 



V 2 



with A$+- and A$,fc given by Eqs. (|79|) and H8U|) respec- 
tively (where, again, care should be taken not to confuse 
/3 as a label on the neutrino flavor with /3 as the source 
angular position). 

Another result of interest is that for the magnification, 
fi^^ , which is defined to be the ratio of the flux of neutri- 
nos of type /3 actually received (from the source at Y and 
given the lensing mass is where it is) to the flux of neu- 
trinos of the same type that would be received with the 
lens absent (but with the source in the same position): 



(89) 



max {{von Blvai B)} = 1. 



This means that 



\N\' 



Npi^thl + 2 J2p,q<p VIpIi 



(84) 



(85) 



Eq. (|83|l establishes the contention made above that in- 
terference effects that emerge with gravitationally-lensed 
neutrinos are a combination of a Young's double slit type 
interference {the [... cos(A$^'')...] envelope term} and flat 
space oscillations [the (Jq,^ — ...) term]. Further, assuming 
the separability requirement is satisfied, we can see how 
\{vi3\ua, Ds)\'^ factorises into an interference pattern and 
a conditional probability. This is a repeat of the behavior 
see in il3.ll 

One should also note that in the particular case of the 
Schwarzschild lens under consideration in the last section, 
the two (assuming non-perfect alignment) classical paths 
from source to detector, denoted by the subscripts -I- and 
experience magnifications given by [see Eq. (2.24) of 

ref. 



(86) 



e xpe 

m 




/++/_+ 2/^ 7^+4^ 



(87) 



where F{Ds, E^) denotes the flux of neutrinos of all types 
that would be received, at an energy of Ei, and factoring 
in geometrical effects, in the absence of the lens. Also, 
Pa denotes the probability that a neutrino generated 
by the source under consideration is of type a. Now, 
for the Schwarzschild lens, assuming mass degeneracy, 
the path difference and phase velocity contributions to 
\{vp\vaTE)s)\'^ factorise into an interference pattern and 
a conditional probability, as mentioned above (Eq. 183|l '). 
This has the effect that the magnification is independent 
of the neutrino flavor under consideration: 



|jVp[/++J_+2^J+J_ cos(A^+-)]P(a^/?)fl,tXP„ 



= [/+ + /_ + 2/r;77cos(A$+-)] 

= ^7^{/5' + 2^|[l + cos(A$+-)]}, (90) 

where P{a — > /3)flat, the flat space neutrino oscillation 
probability, is given by 



P(« ^ /3)fiat 



> ' :sm(A$jfc) 



(91) 



The result for the magnification is as expected given what 
is known about the photon case [see Eq. (9) of Ref. 
Finally for this section, we determine, for future refer- 
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ence, the fringe visibility, V(/3): 
V(/3) 



..max ,,min 



'I'd 



291 



(3^ 



291- 



(92) 



6. PHENOMENOLOGY: HEURISTIC 
CONSIDERATIONS 

Above we have presented the calculation of the phase 
and the consequent phase difference, oscillation probabil- 
ity analog ~ \{vp\vaT Ds)\'^ - which determines the form 
of the oscillation pattern seen at a detector, and mag- 
nification factor, all for the Schwarzschild lens. We now 
turn briefly to the question of the phenomenological con- 
sequences of all these theoretical developments. We shall 
deal with the issues presented here at greater length in 
another work .54] . There are a number of factors which 
broadly determine the visibility of GINI effects [7^ : 

1. Suitability of potentially-lensed sources. The 

first consideration must be, what qualifies as a suit- 
able source? We require sources that produce a 
neutrino signal that might be both gravitationally 
lensed and of sufficient intensity. 

2. Geometrical optics limit. Because our theoreti- 
cal evaluation for the neutrino phase difference has 
been performed within the geometrical optics limit 
(where only the classical paths from source to de- 
tector need be considered in determing the form of 
the interference pattern) , we require that this limit 
holds in the experimental situation under consider- 
ation. This translates to the requirement that the 
de Broglie wavelength of the neutrino mass eigen- 
states in not larger than the Schwarzschild radius 
of the lens (the latter quantity setting the scale of 



the path difference: see below) 



settmg the sc 

ii|n| M- 



3. Detector energy resolution. Even if there ex- 
ists an interference pattern to be mapped out - and 
sufficient events to achieve this - a separate ques- 
tion is whether the smearing of this pattern caused 
by the finite energy resolution of any real neutrino 
detector is so large as to completely wash it out. 

4. Just-so condition for lensing mass. Points [21 
and 131 imply a range for a 'just-right' lensing mass 
(given the energy scale of the neutrinos is already 
set) - not too large and not too small - inside 
which GlNl effects may become evident. This can 
be roughly determined by the following considera- 
tions: for a (point mass) lensing system to produce 
images of similar brightness (so that interference 
effects might be seen), we require that source be 
sufficiently well aligned with the line from the ob- 
server through the lens (i.e., j3 < 9-^). Granted this. 



the scale of the path difference is then set by the 
Schwarzschild radius (see, e.g., [s^l, P- 240), 



Rs = 2GM ~ 3 X 10" 



M 



10-17 Mp 



(93) 



of the lens (and - very broadly - can be consid- 
ered as independent from the distance to the detec- 
tor) , once one has settled on a generic astrophysical 
source which emits neutrinos in some characteristic 
energy range, the lensing mass range is determined. 
This is because we require 



X i?s =i 1 



(94) 



at an energy either within or not too far below 
that detectable by the particular detector technol- 
ogy under consideration (see ii6. 31 below). We label 
this constraint on the lensing mass range the just-so 
condition. 

5. Wave packet spreading and decoherence. By 

analogy with the considerations set out in ^ we 
expect that the full expression for the oscillation 
probability analog include exponential decay fac- 
tors that account for coherence loss effects. These 
essentially factor in the interference attenuation 
which occurs when the different neutrino wavepack- 
ets, traveling with different group velocities and/or 
along paths of different affine length, overlap sig- 
nificantly less than completely at the detector. See 
Hfi.ll for more detail on this issue. 

6. Finite source size effects. Our derivation of 
the phase difference has assumed a stationary point 
source (and detector). Of course, this is at variance 
with Heisenberg uncertainty requirements. But 
more significantly, any real, macroscopic source 
(the region giving birth to all the neutrinos that 
are identified as having come from a particular as- 
trophysical object) will be of finite - indeed macro- 
scopic - size. This can, like detector energy res- 
olution issues, tend to wash away the interference 
pattern because the path difference is now different 
for the various neutrinos that come from different 
parts of the 'same' object. More concisely, an ef- 
fective source angular extent of the order of - or 
larger than - the angular extent of the Einstein ra- 
dius means that the visibility of the interference 
fringes is reduced 0- If the source size is denoted 
by Tsourcej then this translates to the requirement 
that 



< x UGM- 



(95) 



Eq. H95|) is not a sufficient condition to guaran- 
tee a point-like source, however; as energy - and. 
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therefore, phase along any particular path - in- 
creases, there will come a point where (while the 
lens-induced path difference might still generate the 
greatest component of the phase difference for neu- 
trinos from all parts of the lens) the phase differ- 
ence for neutrinos emerging from one part of the 
source will be noticably different to that for neutri- 
nos generated from a different part of the source. 
At this point the interference pattern will, again, 
become smeared out. That this not occur bounds 
the energy: 

E.<^^4^. (96) 

ti source 

7. Finite detector size effects. Much of the dis- 
cussion immediately above carries through, mutatis 
mutandis, to considerations stemming from finite 
detector size. Explicitly, finite detector size effects 
can tend to wash away the interference pattern be- 
cause the path difference (at any particular energy) 
will be non-constant across the volume of the de- 
tector. One must determine whether this is a sig- 
nificant effect. 

8. Finite lens size effects. We have calculated the 
neutrino phase difference in a Schwarzschild metric, 
i.e., assuming the lens to be effectively pointlike. 
This assumption will hold, at least roughly, if the 
Einstein radius of the lensing system [9e, as defined 
in Eq. H73() ] is larger than the scale of the physical 
dimension of the lens. 

9. Source-lens-alignment probability. In order to 
see interference fringes we require that the visibil- 
ity [defined in Eq. H92() ] be sufficiently good. This 
requires a sufficient degree of alignment between 
source, lens, and detector (i.e., a small Y or 77). 
One can then ask, given the lensing mass scale, as 
determined by pointQlabove, and the expected dis- 
tance to a source (of the chosen, generic type), how 
likely is it that there is a lens within a certain dis- 
tance of the line from the source to the detector? 

10. Time scale of lens crossing. Further to the 
point immediately above, one must consider over 
what time scale the lens will cross the 'beam' from 
source to detector and, therefore, how temporally- 
stable - and, indeed, how long-lasting - any inter- 
ference pattern will be. 

11. Intrinsic Source Spectrum. In order to confi- 
dently identify interference effects one must be able 
to rule out the possibility of the intrinsic spectrum 
of the source mimicking these effects. Moreover, 
even given a well-understood source spectrum, a 
separate question is whether there is a measurable 
neutrino flux over a sufficient energy extent that a 
number of interference fringes might be seen at a 
detector. 



6.1. More Detail on Decoherence Effects in Curved 
Spacetime 

We only attempt an heuristic treatment herejsj. Ig- 
noring detector effects (see below), coherence requires 
that there is significant overlap between mass eigenstates 
at a detector. As explained in the appendix, the vari- 
ous mass eigenstates, may have traveled with both dif- 
ferent group velocities and along different paths. Let 
us take a source located on a source plane at Ds and 
neutrino mass eigenstates with an effective width of cTx. 
Then, by analogy with the second exponential damping 
term in Eq. H124|) of the appendix and given the scale 
of the path difference is given by Rs, interference be- 
tween mass eigenstate j traveling down one macroscopic 
path through a Schwarzschild spacetime and k the other, 
roughly requires: 

Here the upper sign refers to the case when the lighter 
mass eigenstate traverses the longer path, the heavier 
along the shorter path, and the lower sign refers to the 
opposite case (there are now four broad cases depending 
on this sign and the relative sizes of Rs and 2^/2ax)- 

Note that if we wish to consider interference between 
different mass eigenstates traveling down the same path 

- i.e., the direct analog of flat space neutrino oscillations 

- we take Rs ^ in in Eq. H97|l so that we require 

< ^^2V2a„ (98) 

(then the equality in the above is satisfied for Ds ~ icoh, 
where Lcoh is the coherence length), whereas if we wish to 
consider interference between the same mass eigenstate 
traveling down different paths, then from Eq. H97(l we 
require 

Rs < 2V2ax. (99) 

Below we shall determine some plausible numbers to put 
in these relations (for a number of different neutrino 
sources), but first we recall some considerations behind 
the determination of ax ■ 



6.2. Determining the Size of the Wavepacket 

At an heuristic level - adequate to the order of magni- 
tude calculations we will make - the neutrino wavepacket 
size (in position space) is given by the size, d, of the re- 
gion to which the neutrino parent particle is localized 

ax ^ d. (100) 
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In turn, d is related to tes, the effective time available for 
the coherent emission, by the parent, of a neutrino wave 
train: 

=^ toff . (101) 

In free space the coherent emission time corresponds to 
the decay time of the parent particle, r, but if the parent 
particle is in a dense and hot medium and undergoing col- 
lisions with its neighbors on a timescale, tcoiiision, smaller 
than T, then tcs ~ tcoiiision |E3j Ell' This effect corre- 
sponds to the collision or pressure broadening of atomic 
spectral lines. In summary, we shall take 

cTx-ltcS with t = min{t coiiiBion,T}, (102) 

where we have explicitly introduced a Lorentz boost, 7, 
to allow for any bulk motion of the source with respect to 
the detector frame. This factor can, of course, be large 
for astrophysical sources. 



are distributed at equal energy intervals, the ab- 
solute uncertainty in neutrino energy determined 
by a detector can be expected to be an increasing 
function of energy. 



7. SUITABLE SOURCE - LENS - DETECTOR 
CONFIGURATIONS FOR GINI 

We can think of four scenarios for source - lens - detec- 
tor configuration that might exhibit GINI effects (there 
may well be more), viz: 

1. Sun - Moon - solar neutrino detector (i.e., in an 
Solar eclipse) 

2. cosmological neutrino source - intervening lensing 
object - large scale Water/Ice Cerenkov neutrino 
detector or airshower array 



6.3. Determing Energy Ranges for GINI 
Phenomenology 

There are two energy ranges that must be considered 
in our analysis, viz 

1. Extrinsic energy range. Forgetting GINI effects 
for the moment, one energy range - which we label 
extrinsic - is delimited by the minimum and maxi- 
mum energies, iJmin and E^ax, at which the generic 
source under consideration can be seen in neutrinos 
by a particular detector technology. These limiting 
energies are determined by either detector or in- 
trinsic source spectrum considerations (whichever 
is the more severe). The extrinsic energy range is 
defined by 

^'min 1^ < ii^max- (103) 



3. artificial neutrino beam on one side of earth aimed 
through center of earth to detector on opposite side 
of the planet 

4. Galactic (i.e.. Milky Way) Core Collapse supernova 
(Types II, lb and Ic) - intervening lensing object - 
solar neutrino detector 

Unfortunately, scenarios 1. to 3. fail one or more of the 
heuristic tests we have set out above and we must, reluc- 
tantly, dismiss them. Scenario 4, however, holds out some 
promise and it is to this that we now briefly turn (see [s^ 
for more detail on all the scenarios mentioned), though 
we alert the reader from the beginning that scenario is 
unlikely to be realised at present because of the low prob- 
ability of supernovae at (neutrino-)detectable distances 
being lensed by objects in suitable mass range. 



2. Intrinsic energy range. We also identify an in- 
trinsic energy range that is given by the follow- 
ing considerations: the lower bound on this range 
is given by the critical energy, -Bcrit, at which the 
pertinent phase difference is equal to one (below 
this value our treatment of the phase breaks down) . 
The scale of this energy is given by requiring 

i^crit ^ 1^. (104) 

Note that the relation is not exact because the RHS 
does not account for the effect of the source align- 
ment parameter, Y , on the phase difference. The 
upper bound on this range, i?washout, is determined 
by the energy at which detector energy resolution 
issues mean that one interference fringe can no 
longer be resolved from another. Washout occurs 
generically because, although interference fringes 



7.1. Core Collapse Supernovae as Sources for GINI 
Observation 

7.1.1. Core Collapse Supernovae as Neutrino Sources: 
Ceneral Considerations 

Let us take the characteristic scale of the distance to a 
Galactic core collapse supernova to be 10 kpc ~ 3 x 10^^ 
cm, the approximate distance to the Galactic Center. 
A core collapse SN observed today at the fiducial 10 
kpc would produce around 10^ and 10'^ events in Su- 
perKamiokande and the Sudbury Neutrino Observatory 
respectively Over the medium term, prospects for 

SN neutrino detection may become even better than at 
present with the construction of the next generation of 1 
Mt underground, water Cerenkov detectors (bOI IgiIIs^. 
For a supernova at 10 kpc, a 1 Mt device should detect 
O(IO^) events 
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7.1.2. Natural Scale for Lensing Mass Required for GIN I 
Effects with SN Neutrinos 

Writing 

A$+- - E,Rs 

we can determine that the smallest lensing mass that 
might produce a phenomenological effect (that we can 
treat using our formahsm) is, very roughly, 10~^^ Mq ~ 
3 X 10^^ g. This is in the cometary mass range. A more 
detailed calculation demonstrates - for the specific 
case of SuperKamiokande - a sensitivity to the GlNl ef- 
fect with lensing masses in the range 

10-1* Mq < Miens < 10'^^ Mq. (106) 

This range is is both conservative and fairly sensitive to 
the SuperK energy thresholds and energy resolution. 

7.1.3. Coherence of Supernova Neutrinos 

A neutrino wavepacket leaving the neutrinosphere of a 
nascent neutron star will have a size [s^ l63j | 

df^-lO-^cm. (107) 

This is to be contrasted with the scale of the affine path 
difference for the lensing mass range under considera- 
tion (Eq. H106|l ^ which is supphed by the range of the 
Schwarzschild radius, viz: 

3xl0~i^cm> Rs >3xl0-"cm. (108) 

We do not, therefore, expect any significant damping of 
the interference amplitude by decoherence due to path 
difference effects. 

There is, however, also decoherence due to group ve- 
locity difference to be considered, i.e., the direct analog 
of decoherence effects for conventional neutrino oscilla- 
tions. The inequality to be satisfied is given by Eq. , 
the RHS of which translates to '--^ 3 x 10^^ cm for 10 
MeV neutrinos jH3 much smaller than the fiducial scale 
of Z?s, ~ 3 X 10^^ cm. We can expect, therefore, to be 
beyond the flat space coherence length. This means that 
the neutrino signal will be characterised by flavor ratios 
that are constant across (measurable) energy. For super- 
nova neutrinos, then, if a suitable lens were present, GINI 
would cause patterns of maxima and minima across en- 
ergy in the detected neutrino spectra. Furthermore, the 
positions, in energy, of these maxima and minima would 
be essentially the same for all neutrino flavors Interfer- 
ence effects would be, in principle, directly evident even 
in neutral current interaction data. On the other hand, 
we would not expect a noticable change in the ratios be- 
tween different neutrino species across energy. In other 



words, for supernova neutrinos, given a suitable lens, it is 
possible to see interference effects due to path difference 
effects but not due to phase velocity difference (i.e., flat 
space oscillation) effects. 



7.1.4. Finite Source Effects with Supernova Neutrinos 

Given a scale for the neutrinosphere, rsN, of ~ 10 
km = 10^ cm, a calculation shows that the point source 
condition, Eq. H95|l . fails at the lower end of the of 
phenomenologically-interesting lensing mass range as- 
suming Dds — Dd- Furthermore, from Eq. (|96|l we find 
that in order that the phase uncertainty introduced by 
the finite size of the supernova neutrino source not be 
too large, we require that the neutrino energy be less 
than ~ 1 MeV, a condition that, with 10 MeV neutrinos, 
we fail to meet by an order of magnitude, again assum- 
ing — D^. We hasten to add, however, that we do 
not believe that either of these two is necessarily fatal: 
a numerical study is needed here and this may well es- 
tablish that GINI effects are visible even when the crude, 
heuristic inequalities above are violated. |83| Moreover, 
that Dds — Dd need not hold (over the Galactic scales 
we are considering) and, further, we might have Ds > 10 
kpc (at the cost of a reduced event rate). In either case 
point source conditions could easily be satisfied. 



7.1.5. Finite Lens Size Effects 

For a lens in the mass range determined above, and 
both source and lens at Galactic length scales, the Ein- 
stein length scale is 

On the other hand, for a lens with the density of ~ 1 gm 
cm-i, the scale of the dimensions of the lensing object, 
I, is given by 

meaning that the classical paths pass very close to the 
object, and, in some cases, one path might actually pass 
inside the object. Given the order of magnitude nature 
of the calculations we have performed here, however, this 
fact will not significantly impinge on the observability of 
the effect we predict. Certainly the neutrinos will not 
interact significantly with the material of the lens. Of 
course, if the lensing object is a black hole, taking the 
lens to be a point source is unproblematic. 
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7.1.6. Finite Detector Size Effects 

A quick calculation shows 0| that finite detector 
size/position resolution effects never become insurmount- 
able over the whole range of possible lens positions. 



1.1.1. Crossing Time Scales for Supernova Neutrinos 

If the lens has a transverse velocity v ~ 30 km s""'^, it 
will move across the Einstein ring in a time scale of ~ 1 
s 0- Given, then, that we expect a detectable neutrino 
signal will be received from a SN over a period of around 
10 seconds, interference fringes will shift over the time 
of observation, but not so quickly that they cannot be 
observed. 



7.1.8. Lensing Probability 

The above paragraphs detail the conditions that a 
source-deflector-observer alignment must satisfy in order 
for GINI to be measured, but implicit at all stages is 
that such an alignment has occured. Unfortunately the 
chance of a suitable deflector lying sufficiently close to 
the line-of-sight to a source in the Local Group is not 
high. 

The Galaxy's rotation curve places a strict upper 
bound on the total mass in its halo (e.g., H^), which 
then implies a maximum possible alignment probability 
to, say, the Magellanic Clouds. Even if the halo consisted 
only of point-masses of suitable size, the simple calcula- 
tion made by Paczynski j65j implies that the lensing opti- 
cal depth, T - essentially equal to the probability that any 
single source is lensed at a given time - of '--^ 10~^. This 
result has been corroborated experimentally by monitor- 
ing stars in both the Galactic center and the Large Mag- 
ellanic Cloud for period variations: both the MACHO 
^ and OGLE jeg groups have found r ~ 3 x lO^^. It 
is important, then, to note that even if the halo is dom- 
inated by point-masses of M ~ 1O~^^M0 suitable for 
GINI, the alignment probability to any neutrino sources 
sufficiently close to be detected at all is only ~ 10~^. 

In the future, however, as detector technology im- 
proves, it may be possible to observe neutrinos from cos- 
mologically distant sources at effective distances of up to 
Gpc. The lensing optical depth is thus increased, both 
because a given mass can act as a more efficient lens and 
because the chance of alignment increases proportionally 
with source distance. A simple calculation of these effects 
(e.g., [s^) imphes that optical depths of close to unity 
are plausible; thus when neutrinos are detected from cos- 
mologically distant sources GINI effects will have to be 
taken into accovmt in the interpretation of any such data 
obtained. 



8. EXTENTION OF THEORETICAL RESULTS 

Besides treating the potential phenomenojogical effects 
of GINI at greater length in another work j5j| , there are, 
of course, a number of directions in which our theoret- 
ical treatment will be extended. Some issues we intend 
dealing with further in another publication include: 

1. From consideration of interference of lensed pho- 
tons in a Schwarzschild metric [see Eq. (7.8) of 
and also see Eq. (9) of 3] Eq. (jSHJ, is actually 
subtly in error: there is an extra —tt/2 phase shift 
missing from the argument of cos term (in other 
words, the interference envelope should actually go 
as -I- 26*1 [1 + sin(A$"' )], generating a central 
minimum for Y — 0). This is present in the case of 
light - and will also be present in the case of neu- 
trinos - because of the opposite parities of the two 
images produced by a Schwarzschild lens (i.e., the 
images - were they able to be distinguished - would 
be flipped with respect to each other: [s^). The 
reason why our treatment has failed to pick this 
extra phase shift up is that we have artificially re- 
stricted the paths under consideration to only the 
classical ones. In other words, we have assumed 
the geometrical optics limit which is strictly only 
valid for phase differences of order unity and larger. 
A more complete treatment using the techniques 
of physical optics - involving integration over all 
paths through the lens plane (each such path be- 
ing uniquely specified by its impact parameter) - 
would recover this phase [and, in fact, demonstrate 
that the full expressions for the oscillation proba- 
bility and magnification involve confluent Hyperge- 
ometric functions that only reduce to trigonometric 
functions in the large phase limit: cf. Eq. (7.11) of 
Ref. [s^l or Eq. (6) of Q]. Moreover, a more com- 
plete treatment would also demonstrate that the 
singularity at /3 = for Eq. (|90|l is not a real effect. 

2. So far we have assumed a static metric. But the 
GINI effect, as noted, requires neutrinos from as- 
trophysical sources that would probably need to 
be located at extra-galactic or even cosmological 
distances for a decent chance that lensing occur 
(though it should be stressed that current detec- 
tor technologies would not allow detection of neu- 
trinos from supernovae beyond out Galaxy and its 
satellites: IsJl) and the introduction of cosmology 
into the formalism developed here would require 
the treatment of a non-static metric £9] . We 
note in passing that were GINI effects ever seen in 
neutrinos from cosmological sources, these effects 
would provide for a test of quantum mechanics over 
the very longest scales. We speculate, then, that 
GINI could be sensitive to the effects of spacetime 
foam (cf., say, |73). In principle, we also expect 
that a GINI pattern in cosmologically-sourced neu- 
trino could be interrogated to determine the value 
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of the Hubble constant Hq 0,|23- Observation (or 
non-observation) of GINI effects would also conti- 
tute a de facto probe of the distribution of dark 
matter objects within a certain well-defined (and 
interesting) mass range jj]. 

3. The Schwarzschild lens is an ideal case never pre- 
cisely encountered in nature. For the realistic situ- 
ation we may need to account for shear multi- 
ple lensing masses, etc. Treating these effects may 
very well demonstrate that the just-so mass range 
is considerably larger than the estimate given by 
Eq. (|106|l 8]. In any case, it is fairly easy to un- 
derstand, at least at the heuristic level, how more 
general lenses might be treated: considering, say, 
Eq. H76|) it can be seen that almost all the con- 
tributions to the phase difference are essentially 
geometrical. The only expression which contains 
information on the mass distribution of the deflec- 
tor is the logarithmic term. This suggests that it 
may be possible to modify our results to arbitrary 
weak deflectors simply by inserting the appropri- 
ate lens potential (and replacing a Schwarzschild- 
specific result for A6'pg). 

4. In we added into the mix coherence loss con- 
siderations. Formally, these can only arise in a 
full wavepacket treatment which we have not at- 
tempted for the curved spacetime case. We remind 
the reader, however, that our results for phase and 
phase-difference will continue to hold in any more 
detailed calculation because these are independent 
of wavepacket considerations. 

5. Our treatment of effects due to the finite nature of 
any real source is very much at an heuristic level. 
Furthermore, finite detector effects can also be im- 
portant, as has been remarked. In this regard, note, 
in passing, that logically connected to this concern 
is the consideration that the observation of GINI 
effects with two - or, preferably, more - widely- 
separated detectors holds out some interesting pos- 
sibilities 0. One would expect here that the in- 
terference patterns seen by different detectors be, 
in general, displaced in energy with respect to each 
other. The degree of this displacement will be re- 
lated to the lensing mass and the geometry. So the 
displacement could probably be used to better con- 
strain relevant parameters than observation with a 
single detector. Moreover, observation of fringes 
with more than one neutrino detector would cer- 
tainly lend credence to the idea that these have 
their origins in GINI. 

9. CONCLUSION 

In this paper we have explicitly calculated the phase for 
a neutrino mass eigenstate propagating through curved 



spacetime, in particular, a Schwarzschild metric. With 
this expression in hand, we have shown how a novel inter- 
ference effect - gravitationally-induced neutrino interfer- 
ence ('GINI') - may show up for gravitationally-lensed, 
astrophysical neutrinos. These interference effects lead to 
a neutrino transition phenomenology qualitatively differ- 
ent from flat space neutrino oscillations. We have shown, 
further, that a result extant in the literature for the 
phase difference with gravitational lensing must be in 
error. We have also derived the form of this phase dif- 
ference when it is given in terms of conventional lensing 
parameters. We have derived the analog of the neutrino 
oscillation probability in flat space for the Schwarzschild 
metric. This quantity controls the phenomenology at a 
detector, in particular, the pattern of maxima and min- 
ima (across energy) for neutrino wavepackets which have 
propagated from source to detector along multiple paths. 
We have adduced heuristic arguments that establish that 
this interference pattern could be seen in the neutrino sig- 
nal from a supernova, provided a suitably-lensed super- 
nova event occurs. Current - and probably even next- 
generation - neutrino detector technologies would seem 
to mean, however, that the probability of such lensing 
occurring for a neutrino-detectable supernova is small. 
Still, for astrophysical neutrinos originating at cosmo- 
logical distances the lensing probability approaches 1 
and some day the technology to detect large numbers 
of these from single sources may become available. We 
have mapped out a program for further research in this 
field. 

In summary, the material presented in this paper serves 
as a proof-of-principle that the GINI effect is both real 
- in a theoretical sense - and, what is more, could lead, 
one day, to interesting phenomenological consequences 
for supernova neutrinos. 
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11. APPENDIX: WAVE PACKET TREATMENT 
OF NEUTRINO BEAM SPLITTER 

In this appendix we set forth a fuU (Gaussian) wave 
packet treatment of the neutrino beam sphter Gedanken 
Experiment (treated in terms of plane waves in §3). Note 
that the results we derive for the exponential damping 
terms (in the equation for the oscillation probability ana- 
log - see Eq. (I124|l ) serve as an heuristic justification of 
the treatment of decoherence we present in §6. 

We write the ket associated with the neutrino flavour 
eigenstate a that has propagated from the source space- 
time position A — [xAitA) to detection position B = 
{xB,tB) as 



p j 

X ldEexp[-i<^P{E;L^^,T^^)]A,{E)\i.,), (Ill) 

where the various quantities are as explained in §3. The 
amplitude for a neutrino created as type a to be detected 
as type P at the spacetime position, B, of the detection 
event is then: 

A, B) = V^E 

p j 

X j dEA,{E)eM~'^'^';{E-L^^,T^^)]. (112) 

Again, we can get rid of the unwanted dependence on 
time by averaging over in the above to determine a 
time-averaged oscillation probability analog at the detec- 
tor position xb ISESES- This gives us that 



\{yp\v^-XA,XB)?^ j dr^^|(i/^|i/a;A,B)p 

pq jk 

X J dEA,{E)AUE)eMKpAE)Lt'''P>^iE)Lf'')], 

(113) 

where one integral over energy has disappeared because 
of the 6{E — E') that arises from the integration over 
time. 

Assuming a Gaussian form for the wavepackets leads 

to 



Aj oc exp 



exp 



4a| 



, (114) 



where Ej is the peak or average energy of mass eigen- 
state J and we employ the very good approximation |44| 
that the wavepacket spread is the same for different mass 



eigenstates. This gives us that 
A,{E)Al{E) 
cx exp 



2£;2 - 2E{E, + Ek) + (-B,' + El) 



4^1 



(115) 



Now defining the peak momentum of mass eigenstate j 
via pj — {Ej ~ mj Y^^, we can write 



where 



p,{E)=p,+Ap,iE), 
Ap,{E) =p,{E)~p, 



(116) 




(117) 



Employing the group velocity of mass eigenstate j, viz. 



■P]/E] 



E^ - m2 ^2 

-~1-^, (118) 



E, 



'1 —3 

we determine that Apj{E) ~ {E — Ej)vj so that 

PjiE) ^ Pj + (E - E,)vj (119) 

We find, then, after a simple calculation that the oscilla- 
tion probability analog becomes 

\{iyf3\l^a;XA,XB)\^ 

- — E v^V? E u^^j^h^pi^^ak 

l^rs ^ pq jk 



X exp 



^{p,Lf^-pkLf^)]} 



X exp 



2 {Ej ~ EkY 



exp 



/„, TAB T AB\2 



(120) 

Note that in the above calculation, though we have been 
employing the group velocity, this does not - and should 
not be seen to - enter into the phase in any fundamen- 
tal way ^3- Indeed, the phase can be calculated en- 
tirely with plane waves (see §3) and, therefore, totally 
without reference to wavepacket notions like group ve- 
locity - though the exponential damping factors above 
critically depend on these. Note also that, again, the 
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normalization has been determined by requiring that 

\{vp\Va\XA,XB)\^ < 1. 

The two damping factors in Eq. Q can be traced back 
to considerations following from (i) source localization 
and the (ii) requirement for overlap of wavepackets at 
the detector's position. Observe that the second damp- 
ing term accounts for an interesting possibility: having a 
heavier - and slower - mass eigenstate travels down the 
shorter path and the lighter mass eigenstate down the 
longer path, will tend to restore coherence. 

As a further particularization of the expression for neu- 
trino oscillation 'probability' with an imaginary beam- 
splitter, we take the expressions for the energies, mo- 
menta, and velocities of the various wavefunctions given 
in terms of expansions around the energy in massless 
limit, Eq^ viz: 



Employing the above we find 



\{iyi3\iya;XA,XB)f 



xexp 



X exp 



X exp 



2 



1 /l^^-L^^ 1 ^]L^^-mlL^^^' 



2ax 



2El 



2a. 



(124) 



E,, 



Pi 



En 



En 



^ ^'2Eo 
^2Eo' 



(121) where the phase difference is given by 



(122) 



where ^ is a dimensionless parameter of order unity de- 
termined by kinematical considerations j40|. The group 
velocity will then be 



4En 



AB 



2En 



(125) 



E, 



2El 



(123) 



Note that we can take the plane wave limit of the above 
equation by setting ^ ^ and Ux ^ oo. This allows us 
to recover Eq. 
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